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SECTION  I 
INTRODUCTION 

The  effects  of  turbulent  media  on  propagation 
studied  for  many  years.  These  studies  have  included  starlight  scintillation, 
underwater  sound  propagation,  laser  beam  propagation  and  satellite  link  scintil- 
lations. Various  theories  and  techniques  have  been  developed  to  estimate  turbu- 
lent effects  under  various  conditions.  The  emphasis  in  this  report  is  on 
satellit e signal  propagation.  Calculational  methods,  particularly  multiple  phase 
screen  algorithms  have  been  developed  that  handle  most  cases  of  interest.  However, 
no  overall  picture  has  been  developed  that  encompasses  the  set  of  likely  in  situ 
environments,  scale  sizes,  and  other  conditions  of  interest.  The  closest  effort 
to  date  was  by  Singleton  (ref.  1),  but  this  work  was  quite  limited  in  the  cases 
studied.  This  lack  has  become  particularly  acute  with  the  growing  volume  of  signal 
data  from  the  Defense  Nuclear  Agency  (DNA)  wideband  satellite  st.jdies.  Air  Force 
Geophysics  Laboratory  (AFGL)  aircraft  studies,  etc.  This  study  is  designed  to 
remedy  the  situation  by  systematically  calculating  and  presenting  in  typical  data 
format  fignal  propagation  results  that  may  be  encountered.  Features  are  identi- 
fied in  the  data  that  allow  interpretation  of  the  in  situ  environments.  Concep- 
tually, the  information  presented  should  be  viewed  as  "clean  data";  i.e., 
experimental  data  without  the  problems  of  lack  of  stationarity,  time-spe'-e  ambigu- 
ities, sampling  problems,  data  sample  boundary  problem:  and  ether  expnritrental 
difficulties  that  becloud  real  data.  One  value  of  point  of  view  s to  realize 
that  features  that  are  marginally  detectable  here  are  probably  not  uso*^^-l  in  real 
data.  Thus,  this  study  defines  the  type  and  qualir.y  of  information  that  can  be 
obtained  from  signal  propagation  data. 

Section  II  discusses  the  propagation  algorithms  pattirularlv  the  Rytov  method 
which  is  used  for  insight  in  interpreting  and  classifying  results.  Section  III 
contains  the  data  analysis.  Section  IV  summarizes  the  conclusions. 

1.  Singleton,  D.  G. , "Saturation  and  Focusing  Effects  in  Radio-Star  and  Satellite 

Scintillations,"  J.  Atoms.  Terr.  Phys.,  Vol  1970. 
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SFCTION  II 

PROPAaATlON  METHODS  AND  SPECTRAL  MODELS 

The  calculations  reported  here  were  done  with  the  multiple  phase  screen  (MPS) 
method  and  the  Rytov  approximation  (RA).  These  methods  are  described  elsewhere 
(refs.  ?.  and  3)  and  only  necessary  extensions  to  these  description  are  given  here. 
The  RA  method  calculates  the  correlation  functions  of  the  log  amplitude  and  phase 
fluctuations  through  a turbulent  medium.  The  basic  two-dimensional  equations  are: 


x(x)  x(x  +'  p) 


(p(x) 


f \z  oos‘[|(z^  * (2) 


(3) 


2.  Wittwer,  L.  A.,  UHF  Propagation  Effects  in  Scintillated  Environments,  AKWL-TR- 
76-304,  Air  Force  Weapons  Laboratory,  Kirtland  AFB,  New  Mexico,  1977. 

3.  Rino,  C.  L.,  Analysis  of  Scintillation  Effects  on  Communication  Systems, 
Stanford  Research  Institute  (to  be  published). 
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where 


z = coordinate  along  propagation  line-of-sight  (i.OS) 

X = coordinate  perpendicular  to  LOS 
Zg  = turbulent  layer  thickness 

Zg  = distance  from  bottom  of  layer  to  receiver  where  the  correlation 
functions  are  measured 
K = rf  wave  number 


Assume  that  there  is  no  variation  of  the  turbulent  layer  in  the  third  coordinate 
and  that  the  fluctuation  statistics  are  isotropic.  The  medium  is  represented  by 


GO  00 

4>2^z,Kpj  = J j An^  (zTx)  An^  (z.e.x  + p)e  ^'^p'^dpde  (4) 

•QO  .00 

where 

An^  * n^  - nj 

n^  » index  of  refraction 
nj  * mean  index  of  refraction 

The  prime  value  of  the  Rytov  equations  is  that  they  provide  a simple,  quanti- 
tatively accurate  estimate  of  the  onset  of  phase  nd  amplitude  fluctuations  and 
their  power  spectra.  If  ^ dependence  then  the  integral  over  z 

can  be  done  analytically,  and  we  nave  for  the  log  amplitude  and  phase  power  spectra 


3 


The  expressions  in  brackets  in  equations  5 through  7 can  be  looked  at  as 
filters  which  filter  the  <I>2  ( j spectrum  to  give  the  different  spectra.  Figure 
1 shews  the  first  two  filter  functions  for  = 10^  cm,  = 3 x 10^  and  K = 0.062. 
The  phase  filter  is  roughly  constant  except  for  some  oscillations.  The  log  ampli- 
tude filter  goes  as  K**  until 

p 


= 5.5  K 
P 


1/2 


(8) 


and  is  approximately  constant  afterward.  This  latter  behavior  can  be  used  to 
evaluate  measured  signal  amplitude  spectra  to  determine  . Before  actual 

data  can  be  confidently  analyzed,  however,  the  practical  limits  of  the  preceeding 
equations  must  be  known,  as  well  as  limitations  on  assuming 


|E(x)|  = 1 + x(x)  (9) 

where  !E(x)|  is  the  field  amplitude.  Equation  9 is  not  actually  necessary,  but 
most  experimentalists  measure  amplitude  and  not  log  amplitude.  The  calculations 
discussed  in  section  III  amply  demonstrate  the  range  of  validity  of  equations  5 
through  7 and  9.  Now  look  at  ^ Assume  that  the  one  dimensional  index 
of  refraction  power  spectral  density  is 


(in) 


where  L is  the  outer  scale  size  and  2v  + 1 is  the  spectral  index.  If  v = 1/2,  we 
have  the  familiar  spectra.  If  the  fluctuations  are  two-dimensional  and  iso- 
tropic, then  the  two-dimensional  power  spectral  density  is 
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Figure  1.  Amplitude  and  Phase  Filter  Functions 
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S^'Vk  K \ = 4tt'/2  ^Ll!ll2L!i^l 

Vp’e)  / \TT1 


f 


l/L*  + K^  + K^ 
P e 


) 


(11) 


For  three-dimensional  isotropic  fluctuations 


(3V  s ^ r (v  + 3/2) 

5i  ) ■ «"  r 

(l/l“  * ♦ K=  * K|  j 


3/2 


(12) 


If  the  fluctuation  statistics  are  not  isotropic,  equations  lo  through  12  can 
be  modified  by  generalizing  the  quadratic  wave  number  terms.  From  equations  (4) 
and  n it  is  clear  that 


Thus  Vj  I Kp j is  the  two-dimensional  index  of  refraction  power  spectral  density  with 
zero  longitudinal  wave  number.  The  one  dimension  correlation  function  correspond- 
ing to  equation  10  is 


An^  (z,'x;  An^  (z.x  + p) 


(14) 


These  formula  are  useable  for  v > - 1/2  but  some  discussion  is  in  order  for 
V < 0.  For  this  case  the  one-dimensional  power  spectral  density  has  infinite 
power,  and  the  correlation  function  is  unbounded  at  p » 0.  The  phase  and  log 
amplitude  spectra,  however,  are  still  finite  because  only  the  spectrum  near  the 
outer  scale  contribute  to  the  propagation.  Thus,  we  can  still  use  an  infinite 
power  spectrum  in  the  sense  that  we  are  fitting  the  real  in  situ  spectrum  correctly 
near  the  sizes  that  determine  the  propagation.  The  use  of  - 1/2  < v _<  0 does  make 
the  total  mean  square  fluctuation  power  ambiguous  because  not  all  of  it  contributes 
to  the  propagation.  Here,  regard  HT  as  just  an  adjustable  parameter 
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K = 0 

” 2r  (}+ Vi)' 


to  determine  by  fitting  to  a measured  power  spectral  density. 
Other  useful  two-dimensional  power  spectral  densities  are 


Chesnut  (ref.  4); 


Gaussian: 


Cylindrical  Rod: 


s5“*(Kp.  * ''c  '• 


|_  V P e / J 

The  one-dimensional  spectra  used  in  this  study  are  shown  in  figure  2.  Other 
useful  functions  for  spectral  modeling  have  been  developed  by  Rino  and  Sachs  (ref. 
5). 

Equations  5 and  6 can  be  used  to  develop  expressions  that  determine  the 
onset  of  amplitude  and  phase  flucturations.  This  was  demonstrated  in  reference  2 
for  the  power  law  expression  and  has  since  been  generalized.  The  following  is 
a summary  for  both  guassian  and  power  law  spectra. 

For 

-r^ 

s 

4.  Chesnut,  W.  G. , "Spatial -Frequency  Analysis  of  Striated  Nuclear  Phenomena— 
Part  2:  A Model  of  the  Striated  Checkmate  Cloud,"  Stanford  Rsch  Institute 
(Unpublished) 

5.  Rino,  C.  L.  and  Sachs,  D.  L.,  "Striation  Models  for  High  Altitude  Nuclear 
Propagation  Effects,"  Stanford  Research  Institute  (to  be  published). 


i 


7 


8 


AFWL-TR-77-183 


For 


<1)2  = x' 


/ 5_.2  An"  L z. 


<l>^  = x' 


/ 2716^ 

\ "><:' , 

1 K'  \ 

L^K 


> 1 


(19) 


(20) 


• (gausslan) 


(21) 


,2  - 


X^  = C(v) 


. / 2Tre^ 

2 

j An*  ir'/*  -i 

(z  + z ) - z* 

vs 

\ mc2 

K"*  L3 

II 

i(f^ 

3 -lv/2  + 1/4 

(power  luw) 


(22) 


(23) 


K 


V+2.5 


V < 1 . 5 


(24a) 


(power  law) 
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x'  = C(v)(^) 


fz  + z 
Is  ( 


V > 1.6 


(24b) 


where 


An  » n - n 

n « electron  density 
n » mean  electron  density 
C(v)  - 0.9  + 0.05 


By  setting  (J>*  ■ x*  * 0.1,  these  expressions  provide  an  accurate  estimate  of  the 
conditions  necessary  for  the  development  of  phase  and  amplitude  fluctuations  for 
uniform  layers. 

Reference  2 demonstrates  that  if  the  scattering  layer  is  sufficiently  thick 
then  the  phase  and  amplitude  distributions  are  a function  of  the  in  aitu  power 
spectral  density.  This  work  is  expanded  to  include  the  phase,  amplitude  and  sig- 
nal power  spectra.  As  in  this  reference,  comparisons  are  made  for  a particular 
in  aitu  spectrum  between  calculations  made  by  assimiing  gaussian  distributed  fluc- 
tuations and  by  using  explicit  geometric  figures. 

The  electron  density  using  explicit  geometric  figures  can  be  written  as 


/ \ ^ 
n(x.z)  ■ E f, 


where 


= random  variable  with  equally  probable  values  of  + 1 
^(^i)  * electron  density 

F ^x,z,x. ,z^ ,a. j = geometric  function  used  to  construct  electron  distribution 
= fluctuation  size 
N = number  of  striations 
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is  used  to  guarantee  that  n(x,2)  is  zero  mean.  The  fluctuation  size,  a^,  is 


distributed  according  to  P(a).  The  P(a)  and  nQ(a)  are  chosen  to  give  the  correct 
in  situ  power  spectral  density  and  An^. 


The  MPS  propagation  requires  the  integrated  phase  of  a plane  wave  going  from 
z to  z + Az  on  a mesh  where  - T < x < T.  In  two-dimensional  calculations 


Z+AZ 


(26) 


For  explicit  representations,  the  striations  are  randomly  located.  If  there 
is  an  average  of  R/A  striations  per  unit  area  in  the  striated  region,  then  each 
phase  screen  will  have  an  average  of  2R  AzT/A  striations  contributing.  The  actual 
number  N for  any  given  screen  is  distributed  according  to  the  Poisson  distribution 


VN) 


.-(i^ 


N 


,.(2N  AzT/A) 


(27) 


To  calculate  a screen,  N is  sampled  along  with  N values  of  a^  and  x^.  The  phase 
screen  is 


2Tre’ 


mc2 


lE'l  "o  F(x.T.x,.a,)<l2 

1 


(28) 


i-0 


The  z limits  are  taken  to  infinity  assuming  that  if  a striation  has  its  center 
within  the  2TAz  area,  then  all  of  the  phase  effects  of  the  striation  are  included 
in  the  phase  screen  being  generated. 


For  gaussian  distributed  fluctuations,  the  phase  screen  is  described  in  terms 
of  a Fourier  series. 


" 2V  ^^1 


-1  iinrx 

n.  ' T- 


(29) 


m=-o 


The  bjjj  coefficients  are  independent,  zero  mean  distributed  variables  with  a 


variance  of 
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= 2T  Az 
m 2 


(?) 


(30) 


where  0^  is  the  two-dimensional  isotropic  in  situ  power  spectrum  defined 
earlier. 


Gaussian,  K*^,  and  cylindrical  rod  power  spectral  densities  are  used  to  illus- 
trate the  convergence  of  the  second  order  propagation  statistics.  The  first  three 
spectra  are  constructed  out  of  gaussian  rods  with 

F(x.j.x,.z,.a,).e'  [(>‘■’'1^  '31) 

The  distribution  functions  for  a and  n^(a)  are 

0 


P(a)  = 6(  L - 2V2  a) 
‘/2 


= /lAni_A  \ 
’ \ N a*  It  / 


gaussian 


(32) 


P(a) 


a > 


21/2  ^2  ^ - 2‘/' 


no(a)  = 


4 An^  A 

NllU*/2  ^ . erf  (V2L)) 


}f^ 


k:' 


..aV4L^ 


(33) 


P{a)  = — t a > 

a'  - 


no(a)  = 


4 An^  A 


ttV2  NJl^  ( 1 - erf  (v2l)) 


1/2 


K-* 


..aV4L2 


J 


(34) 


If, 

“j 

\ 

♦} 
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where  «.  is  the  inner  scale  size,  6 (a)  is  the  Dirac  delta  function,  and 


erf  (x)  = ^7T  J e"^  dt 
0 


(35) 


The  cylindrical  rod  function  is 


’ V 2 

/ \ ^ 

= 1 1 

X - Xi 

t(z-Zi) 

L ' 

nl/2 


< L 


(36) 


and 


cylindrical  rod 


(37) 


The  measure 


P{a)  » 6 (l  - 2'/^  a) 

V J/2 

, s { V 

n„(a)  * I — r — ) 

° \ Zv  n f 

of  thickness  of  a scauering  layer  is  the  mean  free  path  defined 


as 


MFP 


J (2a)  f P(a)  da 


t/2*/‘ 


(38) 


for  distributions  where  P(a)  ^ 6 (l  - 2‘/^  a)  and 


MF? 


[^] 


(39) 


for  the  others. 

For  a given  spectral  type.  z^.  and  L.  N/A  is  varied  to  give  the  number 
of  mean  free  paths  desired.  The  calculations  described  in  section  HI  examine 
convergence  for  the  second  order  signal  statistics  with  scattering  layers  that 
are  five  mean  free  paths  thick. 
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SECTION  III 
DATA  ANALYSIS 

Calculations  were  run  for  five  in  situ  spectral  types;  K"*,  K"*,  K"^  K”®, 
and  gaussian.  Data  taken  by  AFGL"  and  SRI**  (the  latter  from  the  DNA  Wideband 
Satellite)  imply  the  existence  of  a K"*  spectrum  in  polar  environments  .with  a 
reasonable  frequency  of  occurrence.  Goldman  (ref.  6)  has  derived  conditions  under 
which  this  spectrum  might  be  expected.  The  spectrum  has  been  observed  in 
several  ionospheric  experiments  (refs.  7 through  9)t.  This  spectrum  is  considered 
the  result  of  the  IxF  instability  (refs.  6,  10  through  12).  Recently  it  has  been 
shown  that  K“*  can  result  from  the  decay  of  FxF  modes  into  drift  modes  (ref.  13). 
The  K"®  spectrum  also  has  been  inferred  from  AF6L  and  SRI  data.  This  spectrum 
has  been  calculated  at  NRL  (ref.  14)  in  spread  F simulation  studies.  The  K~®  and 
gaussian  spectra  were  used  to  examine  the  propagation  effects  of  single  size 
spectra,  which  have  not  appeared  commonly  in  data. 

In  the  following,  the  environments  are  described  by  specifying  L.  z..  z .K. 

— ™ 5 g 

An*  and  the  type  of  fluctuation  power  spectral  density.  The  propagation  LOS  is 
assumed  to  travel  through  the  stationary  fluctuations  at  v » 3 x lO**  cm/sec. 

This  velocity  is  used  to  transform  spatial  data  to  time  data.  In  particular,  to 
analyze  the  calculated  spectral  data,  it  is  convienent  to  define  three  frequencies 
which  correspond  to  the  first  Freznel  zone  size,  the  flucturation  scale  size  and 
the  signal  scale  size,  respectively. 


* Whitney,  H.,  Air  Force  Geophysics  Laboratory,  private  communication. 

**  Rino,  C.  L., Stanford  Research  Institute,  private  communication. 

t Kelley,  M.  C.,  Cornell  University,  private  communication. 

NOTE:  References  6 through  14  are  listed  at  the  back  of  this  section. 
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wry 

(power  law) 

(41a) 

"my 

(gaussian) 

(41b) 

1 

7 

(42) 

po) 

The  Fresnel  frequency,  fp,  corresponds  approximately  to  the  maximum  in  the  ampli- 
tude filter  function  transformed  to  frequency  space. 

There  are  two  classes  of  fluctuation  power  spectra  as  distinguished  by  the 
signal  structures  they  give  rise  to.  The  first, called  multisize  spectra, includes 
all  the  power  law  spectra  for  v < 1.5.  The  name  derives  from  the  observation 
that  the  size  of  fluctuation  that  dominates  the  amplitude  fluctuations  for  weak 
scatter  and  fj  1 fp  is  determined  by  the  Fresnel  zone  size  and  not  the  fluctuation 
scale  size.  Thus,  any  size  fluctuation  can  dominate  the  amplitude  effects  given 
the  proper  circumstances.  The  second  class,  called  single  size  spectra,  incluoes 
power  law  with  v < 1.5,  gaussian,  and  exponential  power  spectral  densities.  For 
these  cases  the  amplitude  effects  are  dominated  by  a single  size  fluctuation,  the 
scale  size  L. 

For  future  purposes,  it  is  useful  to  distingui^  between  weak  and  strong 
scatter  cases.  Weaj^scatter  is  defined  as  having  1 0.1  and  ^ 0.5  while 
strong  scatter  has  > 0.1  or  ^ 0.5.  The  definition  for  weak  scatter  has 
been  chosen  because  as  demonstrated  by  the  calculations,  the  weak  scatter  region 
is  where  equations  5 through  7 and  9 are  applicable. 

The  interplay  of  environmental  parameters  defines  different  types  of  signal 

properties  for  each  spectrum.  This  is  demonstrated  by  the  propagation  sp  ce 

plots  in  figures  4,  20,  39,  55  and  72  which  were  constructed  using  equations 

19  through  24  to  define  x^  = 0.1  and  4)^  = 0.1  curve".  In  all  but  the  gaussian 

plot,  z ,z  ,K  and  spectrum  type  are  fixed,  where  An^  and  L are  the  varied  quanti- 
5 g — 

ties.  The  gaussian  case  fixed  An  but  varied  z^  with  all  other  variables  the  same 
as  the  other  spectra  calculations  to  examine  the  z^  dependence  in  equation  23. 

Each  of  the  propagation  space  plots  has  .three  distinct  regions  shown  schematically 
in  figure  3. 
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L(cm) 

Figure  3.  Propagation  Space  Plot  fo;  ''"neric  Spectrum 

a.  Region  I is  characterized  by  gaussian  quadrature  evolving  to  Raleigh 
statistics  as  An^  or  increases  and  by  fp  f^.  The  sets  of  figures  7 through  9, 
42  through  44,  59  through  61  and  75  through  85  show  the  evolution  of  the  signal 
statistics.  It  is  evident  from  the  complex  signal  plots,  the  amplitude  distribu- 
tions, and  the  phase  distributions  that  the  statistics  are  approximately  gaussian 
quadrature.  It  is  also  clear  that  the  distribution  plots  are  not  very  sensitive 
to  in  situ  spectral  type.  The  different  spectra  do  manifest  themselves  in  the 
complex  signal  plots  and  the  spectral  densities.  Figures  23a  and  24a  for  the 
K"^  spectrum  and  figures  75a  and  76a  for  the  gaussian  spectrum  clearly  show  the 
relative  lack  of  high  frequency  fluctuation  in  the  gaussian  case.  The  amplitude 
and  phase  spectra  for  the  weak  scatter  cases  are  clearly  explained  by  the  Rytov 
approximation  equations.  Figures  23  and  75  are  representative  of  fp  < f^  where 
the  peak  in  the  amplitude  filter  function  occurs  before  the  break  in  at  f^. 
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The  oscillatory  behavior  of  the  filter  functions  is  apparent.  The  amplitude 
spectra  begin  with  an  approximate  f**  behavior  to  fp  and  the  behavior  after  fp. 
For  power  law  spectra  with  a characteristic  well  defined  break  at  f^,  the  region 
between  f^  and  fp  will  be  approximately  flat.  The  phase  spectrum  for  the  most 
part  mirrors  $2  with  the  phase  filter  function  effects  apparent.  The  Fresnel 
frequency  is  denoted  by  a deep  notch.  In  some  cases,  the  notch  is  easier  to 
identify  than  the  corresponding  amplitude  spectrum  peak. 

If  fp  can  be  found  by  examining  either  the  amplitude  or  phase  spectrum 
and  it  is  at  the  low  frequency  end  of  the  flat  portion  of  the  amplitude  spectrum, 
then  fp,  fg  and  the  spectral  type  are  unambiguously  determined.  The  location  of 
fp  must  be  identified  because  the  same  spectrum  behavior  can  be  generated  with  a 
power  law  spectrum  with  v » 1 where  > ^5*  and  fp  will  be  at  the  high  frequency 
end  of  the  flat  portion  of  the  amplitude  spectrum. 

Even  if  fp  is  not  obvious  from  the  spectrum,  often  a good  estimate  can 
be  obtained  from  knowledge  of  the  propagation  geometry  and  interpretation  anibigu> 
ities  can  be  reduced.  Figures  7,  26,  42,  59  and  81  show  several  weak  scatter 
spectra  where  f^  » fp.  The  flat  portion  of  the  amplitude  spectrum  has  disappeared 
and  the  spectrum  now  has  a sharp  peak,  let  us  define  the  difference  between  tie 
approximate  power  law  exponent  of  the  amplitude  spectrum  for  f < f^  as  o^. 

For  the  spectra  considered  here,  ^ 5 for  these  peaked  spectra. 

This  behavior  unambiguously  defines  f^,  fp  and  the  spectral  type.  These  sharp 
peaked  data  have  been  measured  by  AFGL.* 

For  all  of  the  weak  scatter  cases  with  fp  f^,  the  signal  power  spectrum 

is  approximately  the  sum  of  the  amplitude  and  phase  spectra.  This  is  seen  by 
examining  equations  5 through  7 and  noting  that  ^ » E|  and  « Ej- 

As  An^  or  increases,  the  scattering  becomes  strong.  The  phase  spectrum 
becomes  ambiguous  because  of  the  rapid  phase  changes  during  deep  fades.  The  anqili' 
tude  spectrum  now  begins  to  distort  in  shape  where,  for  weak  scatter,  it  just 
increased  in  magnitude.  The  changes  start  at  the  lowest  frequency  and  then  work 
upward.  The  general  trend  is  toward  a flat  spectrum  with  a break  at  f^  and  with 
the  $2  behavior  beyond.  Figure  27d  shows  a phenomenon  that  we  will  often  see  in 


* Whitney,  H.,  Air  Force  Geophysics  Laboratory,  private  communication. 
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in  strong  scatter  for  all  types  of  propagation.  There  is  an  apparent  second 
break  in  the  amplitude  spectrum  that  does  not  correspond  to  any  characteristic 
frequency.  This  break  designates  the  failure  point  in  frequency  below  which  the 
Rytov  equations  no  longer  hold.  Thus,  equations  fail  first  ac  small  f,  and  the 
point  of  failure  moves  to  higher  frequency  in  a well-defined  way  eventually 
coinciding  reasonably  well  with  f^.  Behavior  such  as  in  figure  27d  has  been 
observed  in  data  taken  by  AFGL.* 

The  complex  signal  spectrum  does  not  change  shape  appreciably  as  An^ 
increases  but  the  break  at  f^  does  move  to  higher  frequency.  The  spectral  shape 
after  f^  tends  to  follow  tj.  Figures  6,  22,  41,  58  and  74  show  plots  of  versus 
X*  for  all  of  the  spectra.  The  is  approximately  the  same  monotonic  function 
of  X*  regardless  of  in  aitu  spectrum  when  f > i, , 

5 C 

b.  Region  II  is  characterized  by  doughnut  statistics  evolving  to  Rayleigh 
as  An^  increases  and  by  f^  > fp.  Doughnut  statistics  are  characterized  by  large 
phase  fluctuations  and  small  amplitude  fluctuations.  These  statistics  are  often 
seen  In  the  ambient  ionosphere  and  are  thought  to  be  characteristic  of  other 
disturbed  ionospheres  as  well.  The  general  behavior  in  region  II  is  a function 
of  whether  the  in  aitu  spectrum  is  multisize  or  single. 

(1)  Multi  size  Spectra.  Figures  30,  31,  32  and  33  show  a typical 
sequence  of  multisize  calculations  from  the  K“*  data.  The  first  two  sets  of 
figures  are  weak  scatter  and  demonstrate  the  preponderance  of  phase  effects. 

The  dominance  of  phase  effects  increases  with  L as  demonstrated  by  comparing 
the  above  figures  with  figures  34,  35  and  36.  The  influence  of  the  particular 
type  of  multisize  in  aitu  spectrum  is  apparent  in  the  complex  signal  plots  seen 
in  figures  31  and  11  which  correspond  to  environments  that  are  significantly 
different  only  In  the  in  aitu  spectrum  assumed. 

Figures  35d  and  35e  show  good  weak  scatter  examples  of  phase  and 
amplitude  for  multisize  spectra.  The  phase  spectrum  closely  maps  except  for 
the  notch  at  fp  which  represents  the  power  lost  to  the  amplitude  fluctuations. 

The  amplitude  spectrum  demonstrates  the  behavior  expected  from  examining  the 
Rytov  equations.  There  are  three  distinct  spectral  regions.  The  first  is  for 
f < f^  and  is  driven  by  the  f behavior  of  the  amplitude  filter  function  because 

* Whitney,  H.,  Air  Force  Geophysics  Laboratory,  private  communication. 
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$2  is  approximately  flat.  For  ^5  < f < fp»  'I’j  “ the  filter  function  still 

goes  as  f resulting  in  an  approximate  f‘  behavior.  For  f > fp,  the  filter  func- 
tion becomes  flat  and  the  amplitude  spectrum  goes  as  $2*  the  difference,  a^,  in 
the  power  law  exponent  for  f > f^  and  f < f^  is  approximately  4 regardless  of 
the  spectral  type.  It  is  clear  that  the  amplitude  spectral  shape  near  the  peak 
is  only  a function  of  the  spectral  type  and  not  L. 

Since  the  total  amplitude  power  is  dominated  by  the  spectrum  near 
the  peak,  amplitude  for  multisize  spectra  can  be  characterized  by  the  spectral 
type  and  x*-  This  is  an  important  result.  It  says  that  if  lines  are  drawn  in 
region  II  parallel  with  the  ® 0.1  curve,  the  amplitude  spectrum,  power,  and 
distribution  are  constant  along  any  of  the  lines.  This  is  very  useful  in  that 
a calculation  for  a given  x^  defines  the  amplitude  effects  for  all  L and  An* 
with  the  same  x^  with  Tp  < f^.  This  means  that  from  amplitude  data  alone,  only 
the  distance  away  from  the  ^ ^ 0.1  line  is  determinable  and  not  the  position 
along  that  line.  Thus,  only  An*  L"*'^  can  be  determined.  If  L or  An*  can  be 
deduced  from  the  phase  or  other  data,  then  the  ambiguity  can  be  resolved.  The 
complex  signal  spectrum  for  ^ < 0.1  is  prooortional  to  As  the  propagation 
effects  increase,  the  break  denoted  by  f^^  moves  to  higher  frequency  keeping  the 
«'2  behavior  for  f > f^.  These  calculations  show  that  both  and  fp  are  usually 
apparent  in  the  weak  scatter  amplitude  and  phase  spectra  as  long  as  the  frequency 
range  is  large  enough.  Past  scintillation  data  have  tended  to  emphasize  resolu- 
tion at  ever  higher  frequencies.  The  low  frequency  data  are  also  quite  valuable 
particularly  to  find  f^  even  when  the  corresponding  wavelengths  may  be  as  large 
as  the  characteristic  dimension  of  the  ionosphere.  As  An*  increases,  x^  becomes 
larger  than  0.1,  becomes  larger  than  0.5,  and  the  propagation  signal  scatter- 
ing becomes  strong.  The  phase  spectrum  for  all  in  situ  spectra  becomes  dominated 
by  large  rapid  phase  changes  and  no  longer  accurately  maps  4^^.  The  amplitude 
spectrum  evolves  qualitatively  as  in  region  I as  the  scattering  becomes  strong. 

The  Rytov  equations  begin  to  fail  starting  at  the  lower  frequencies.  The  tendency 
is  for  the  low  frequencies  to  flatten  out  to  a break  at  approximately  f^  and 
behavior  beyond  f^.  See  figures  JOd,  31d,  32d,  and  33d.  The  amplitude  statistics 
continue  to  be  constant  along  lines  parallel  to  x‘^  =0.1  in  the  strong  scatter 
region.  The  complex  signal  spectrum  changes  are  characterized  by  the  increases 
ifi  with  little  change  in  spectrum.  The  and  x^  increase  essentially  mono- 
torically  to  limits  of  1.0  and  0.46,  respectively.  The  amplitude  distributions 
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evolve  to  Rayleigh  showing  little  spectral  dependence.  In  addition,  the  ampli- 
tude distributions  always  have  a smaller  population  of  deep  fades  than  a Rayleigh 
distribution  (figures  37b,  32b  and  53b).  This  is  the  source  of  the  statement 
that  Rayleigh  statistics  are  a worst  case  for  amplitude  sensitive  systems.  The 
discussion  of  single  size  spectra  shows  that  this  is  not  always  true. 

For  both  weak  and  strong  scatter  is  the  same  well  defined,  single 
valued  function  of  that  was  found  in  region  I.  From  the  discussion  of  the 
amplitude  spectrum,  it  is  apparent  why  this  is  true  for  any  given  spectral  type, 
but  it  is  not  obvious  why  the  same  function  works  for  all  multisize  spectra  and 
for  all  spectra  in  region  I.  The  most  probable  cause  is  that  (in  the  cases 
quoted)  the  majority  of  amplitude  power  comes  from  a narrow  spectral  range. 

(2)  Single  Size  Spectra 

The  behavior  of  single  size  spectra  in  region  II  is  significantly 
different  and  distinctive  from  that  of  the  multisize  spectra.  The  complex  signal 
plots  clearly  lack  the  high  frequency  content  of  the  multisize  spectra  (figures 
64a,  65a,  84a  and  85a).  For  weak  scatter,  the  amplitude  and  phase  spectra  are 
precisely  what  would  be  expected  from  the  Rytov  equations.  The  basic  difference 
from  earlier  cases  is  that  now  breaks  faster  than  f*  and  the  amplitude  spec- 
trum peak  is  now  determined  by  the  scale  size  and  not  the  Fresnel  zone  size. 

Figures  86  through  90  show  a typical  sequence  of  calculations  as  An^  increases 
for  the  gaussian  spectrum.  Figures  68  through  71  show  weak  and  strong  scatter 
cases  for  the  spectrum.  These  figures  clearly  demonstrate  the  folding  of  the 
filter  functions  with  <l>2.  In  figures  68d,  the  amplitude  spectrum  goes  as  f*  for 
f < r,  f”*  for  f-  < f < fc.  and  f"®  for  f > f^.  The  Rytov  equations  show  that 
for  any  v > 1.5  there  is  some  f^  < fp  for  which  that  portion  of  the  amplitude 
speitrwn  containing  most  of  the  power  goes  as  K"**  (k^)  independent  of  fp. 

The  size  of  the  ratio  fp/fj  required  to  make  the  amplitude  spectrum  and,  hence, 
the  amplitude  statistics  independent  of  fp  depends  upon  how  large  v is.  In  the 
and  the  gaussian  calculations,  we  just  barely  reach  the  region  of  fp  indepen- 
dence as  judged  by  the  amplitude  spectrum  before  the  numerical  limits  of  the 
propagation  algorithm  are  reached.  For  very  small  f^  compared  to  fp  then,  we 
see  that  the  amplitude  statistics  and  spectrum  are  constant  on  curves  parallel 
to  the  =0.1  curve.  We  also  see  that  the  versus  curve  becomes  invariant 
along  = constant  and  roughly  independent  of  single  size  spectral  type.  This 
function  is  different  from  the  corresponding  function  for  multisize  spectra  and 
offers  an  easy  way  to  partially  classify  spectra  from  and  data  (figures  * 

I 
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58  and  74).  It  must  be  remembered,  however,  that  must  be  sufficiently  small 
compared  to  fp  for  the  two  types  of  spectra  to  be  distinguishable  in  this  manner. 

Figures  56  and  72  show  the  propagation  space  plots  for  the  r®  and 
gaussian  spectra,  respectively,  with  = 0.5  plotted  in  addition  to  the  = 0*1 

curve.  The  criterion  for  weak  scatter  is  now  driven  by  alone.  The  larger 

for  a given  compared  to  a multi  size  case  is  evident  from  the  complex  signal 
plots.  These  plots  are  characterized  by  a lack  of  deep  fades  and  large  signal 
enhancements  (figures  65a,  70a,  71a,  88a,  92a,  65b,  70b,  88b  and  92b).  It  is 
these  enhancements  that  drive  preferentially  to  x^-  The  amplitude  distribu- 
tions also  illustrate  the  presence  of  large  enhancements  and  low  deep  fade  popu- 
lations. 

As  An^  or  z^  increases  and  the  scattering  goes  from  weak  to  strong, 
the  spectra  no  longer  obey  the  Rytov  equations.  The  phase  spectrum  becomes  domi- 
nated by  passes  near  the  origin  and  becomes  useless.  The  amplitude  spectrum  tend 
to  distort  and  change  at  all  frequencies,  not  just  the  low  frequencies  as  in  the 
multisize  spectra.  A new  break  appears  which  is  correlated  with  f^  (figures  66d, 

71  d,  84d,  87d,  89d  and  92d).  The  slope  for  ^ > i^^  tends  to  a 5.0  to  5.5  power 

law  for  both  the  gaussian  and  the  K’®  spectra. 

The  intermediate  slope  has  a range  of  slopes.  The  most  interesting 
feature  of  the  strong  scatter  is  that  becomes  significantly  greater  than  one 
driven  by  the  large  signal  enhancements.  The  peak  in  the  amplitude  distribution 
moves  to  deeper  fades  and  soon  this  case  poses  a more  severe  fade  threat  than 
Rayleigh  in  having  more  deep  fades.  Eventually  as  An^  increases,  peaks,  then 
goes  to  1;  and  the  amplitude  distribution  goes  to  Rayleigh.  The  complex  signal 
spectrum  shows  a unique  behavior  not  seen  in  multiple  size  spectra.  At  large 
frequencies,  the  behavior  maintains  itself.  At  the  lower  frequencies,  however, 
the  spectrum  is  evolving  to  a gaussian.  Figures  69f  and  66f  demonstrate  this 
behavior.  Eventually,  the  complex  signal  spectrum  will  be  gaussian  for  all 
important  frequencies. 

Figure  33a  shows  a complex  signal  plot  for  a case  that  has  Rayleigh 
statistics;  figure  85a  is  a corresponding  Rayleigh  case  for  the  gaussian  spectrum. 
Even  though  both  are  Rayleigh,  the  difference  in  high  frequency  content  is  evident. 
This  is  traceable  to  a K'^  complex  signal  spectrum  in  one  case  and  a gaussian  in 
the  ocher.  It  seems  clear  that  eyeball  analysis  of  signal  plots  is  a useful 
technique  for  differentiating  types  of  spectra. 
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c.  Region  III  is  the  limiting  statistics  of  the  assumed  spectra.  The 
amplitude  statistics  are  Rayleigh.  The  amplitude  spectrum  is  flat  at  f < f^  and 
drops  off  rapidly  above  f^.  The  complex  signal  spectrum  evolves  to  an  approxi- 
mate gaussian  for  all  cases  except  for  where  the  asymptotic  strony  scatter  limit 
is  approximately  K"^.  The  shapes  of  these  spectra,  except  perhaps  for  the  K"^, 
do  not  contain  much  useable  information  about  the  environment.  The  only  para- 
meter remaining  which  characterizes  both  the  complex  signal  and  amplitude  spec- 
tra is  Jlj.  Figures  5,  21,  40,  57  and  73  shew  plots  of  for  all  the  studied 
spectra.  The  Lj^  and  1,^  are  the  environment  scale  sizes  perpendicular  and  paral- 
lel, respectively,  to  the  propagation  LOS.  The  plots  are  similar  except  for  the 

K"*  case. 

p 

Figures  93  through  96  show  the  results  comparing  the  signal  statistics 
for  gaussian  distributed  fluctuations  with  those  from  an  explicit  representation 
of  the  striations  for  four  different  in  situ  power  spectral  densities.  The 
explicit  representation  results  are  represented  by  dashed  histograms  on  the  dis- 
tribution plots  and  circles  on  the  power  spectrum  plots.  Table  6 shows  various 
other  statistical  measures  for  comparison.  All  of  the  calculations  with  explicit 
geometric  representations  are  five  mean  free  paths  thick. 

Figures  93  and  94  represent  the  gaussian  and  cylindrical  rod  in  situ 
power  spectra,  respectively.  The  agreement  is  evident  between  the  gaussian  dis- 
tributed fluctuations  and  the  explicit  representations.  Figures  95  and  96  are 
for  the  Kp*  and  K"‘  cases,  respectively.  Again,  the  agreement  is  excellent  up 
to  about  3 hertz.  The  lack  of  agreement  above  3 hertz  results  from  the  minimum 
size,  £,  used  in  the  explicit  representations,  because  l was  set  to  2 x 10^  cm 
which  corresponds  roughly  to  3 hertz.  Table  6 shows  that  the  agreement  extends 
to  other  important  statistical  measures.  These  results  demonstrate  conclusively 
that  for  layers  that  are  five  or  more  mean  free  paths  thick  the  in  situ  power 
spectral  density  is  necessary  and  sufficient  for  calculating  the  required  signal 
statistics.  The  assumption  to  gaussian  fluctuations  for  the  use  of  explicit 
representation  is  merely  a matter  of  convenience. 
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SECTION  IV 
CONCLUSIONS 

This  report  illuminates  the  basic  character  of  propagation  through  turbulent 
environments  characterized  by  different  in  aitu  power  spectral  densities.  The 
Rytov  method  was  used  to  develop  analytic  comprehendabl e expressions  for  the 
phase  and  amplitude  power  spectral  densities,  equations  5 and  6,  for  relatively 
uniform  layers.  Expressions  (equations  19  through  24)  were  also  developed  for 
the  mean  square  log  amplitude  and  phase  fluctuation  power.  Calculations  using 
the  MRS  algorithm  shows  that  these  expressions  hold  quite  well  as  long  as 
^ < 0.1  and  < 0.5,  subsequently  defined  as  weak  scatter  criteria,  even  witr. 
the  additional  assumption  that  |E(x,z)|  » 1 -«■  x (x.z).  Numerous  examples  show 
that  features  exist  in  weak  scatter  spectra  that,  via  the  Rytov  expressions, 
allow  determination  of  scale  sizes,  Fresnel  zone  sizes,  spectral  indices,  and 
other  valuable  information. 


For  any  given  in  aitu  power  spectral  density,  three  classes  of  propagation 
effects  exist  which  correspond  to  the  three  regions  of  a propagation  space  plot. 
The  distinctive  character  of  each  of  these  regions  is  readily  apparent  in  the 
signal  statistics.  In  region  II,  distinguished  by  f^  < fp,  the  study  shows  that 
the  general  behavior  of  the  propagation  is  a function  of  whether  the  in  aitu 
spectrum  is  single  or  multisized.  The  multisize  effects  are  dominated  by  the 
first  Fresnel  zone  size,  while  the  single  size  cases  are  dominated  by  the  scale 
size.  In  both  cases  for  fp  sufficiently  large  compared  to  f^,  the  amplitude 
statistics  become  constant  along  lines  parallel  to  the  = 0.1  curve  in  the 
propagation  space  plots  for  both  weak  and  strong  scatter.  This  result  is  very 
important  in  that  a single  set  of  calculations  can  characterize  the  amplitude 
effects  for  all  f^  < fp.  The  two  types  of  spectra  are  further  differentiated 
by  their  plots  of  versus  x^*  Th®  differences  in  these  curves  are  sufficient 
to  constitute  a very  good  test  as  to  spectrum  type.  In  addition,  the  curve  for 
the  multi  size  spectrum  holds  for  all  spectra  for  f^  ^ fp. 

All  of  these  results  imply  that  for  the  best  environmental  analysis,  weak 
scatter  cases  are  best  at  the  highest  propagation  frequencies  possible  to  maximize 
the  fluctuation  power  being  analyzed.  A number  of  calculations  examine  how  thick 
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a turbulent  layer  has  to  be  so  that  the  propagation  effects  only  depend  uniquely 
on  the  in  situ  power  spectral  density  and  the  propagation  geometry.  Results  show 
that  five  mean  free  paths  are  sufficient  to  minimize  any  dissimilarities  between 
different  electron  density  representations  which  have  common  power  spectral 
densities. 


24 


AFWL-TR-77-183 


REFERENCES 

1.  Singleton,  D.  G. , "Saturation  and  Focusing  Effects  in  Radio-Star  and 
Satellite  Scintillations,"  J.  Atmos.  Terr.  Phys.,  Vol  1970. 

2.  Wittwer,  L.  A.,  UHF  Propagation  Effects  in  Scintillated  Environments,  AFWL-TR- 
76-304,  Air  Force  Weapons  Laboratory,  Kirtland  AFB,  NM,  August  1977. 

3 . Ri no , C . L . , Analysis  of  Scintillation  Effects  on  Communication  Systems, 
Stanford  Research  Institute  (to  be  published). 

4.  Chesnut,  W.  G. , "Spatial -Frequency  Analysis  of  Striated  Nuclear  Phenomena-- 
Part  2:  A Model  of  the  Striated  Checkmate  Cloud,"  Stanford  Research  Institute 
(unpublished). 

5.  Rino,  C.  L.  and  Sachs,  0.  L.,  "Striation  Models  for  High  Altitude  Nuclear 
Propagation  Effects,"  Stanford  Research  Institute  (to  be  published). 

7.  Phelps,  A.  D.  R.  and  Sagalyn,  R.  C.,  "Plasma  Density  Irregularities  in  the 
High-Latitude  Top  Side  Ionosphere,"  of  Geophysical  Rsch,  Vol  No.  4, 
February  1976. 

8.  Dyson,  P.  L.,  McClure,  J.  P.,  and  Hanson  W.  B.,  In  Situ  Measurements  of  the 
Spectral  Characteristics  of  F Region  Ionospheric  Irregularities,  Vol  79, 

April,  1974. 

9.  Morse,  F.  A.,  et  al.,  EQUION,  An  Equatorial  Ionospheric  Irregularity  Experi- 
ment,  SSL-76  (6960-04)“2,  Aerospace  Corporation,  El  Segundo,  Calif.,  April 
1976. 

10.  Scannapieco,  A.  J.,  Ossakow,  S.  L.,  Goldman,  S.  R. , and  Pierre,  J.  M. , 

"Plasma  Cloud  Late-Time  Striation  Spectra,"  J.  Geophys.  Rsch. , Vol  1976. 

11.  Ragnlien,  T.  0.  and  Weinstock,  J.,  "Theory  of  the  Nonlinear  Spectrum  of  the 
Gradient  Drift  Instability  in  the  Equatrial  Electrojet,"  J.  Geophys.  Rsch. . 

Vol  79,  November  1974. 

12.  Sleeper,  A.  M.  and  Weinstock,  J.,  "Nonlinear  Theory  of  Density  Fluctuation  in 
Turbulent  Plasmas,"  Phys.  Fluids,  Vol  1^,  1972. 

13.  Chaturvedi,  P.  K.  and  Kaw,  P.  K. , "An  Interpretation  for  the  Power  Law  Spectrum 
of  Spread  F Irregularities,"  J.  Geophys.  Rsch.,  Vol  1976. 

14.  Scannapieco,  A.  J.,  Goldman,  S.  R.,  Ossakow,  S.  L. , Bock,  D.  L.  and  McDonald, 

B.  E.,  "Threoretical  and  Numerical  Simulation  Studies  of  Midlatitude  F Region 
Irregularities,"  NRL  Memorandum  Report  3014,  Naval  Research  Laboratory, 

Wash  D.  C. , March  1975. 


25/26 


APPENDIX 

CALCULATION  DATA 

This  ai^pendix  contains  the  sets  of  data  which  different  in  situ  power  spectra, 
L,  and  propagation  geometry.  All  of  the  striated  layers  are  uniform  in  their 
statistical  properties. 

The  signal  incident  on  a striated  layer  beginning  at  z = 0 is  a plane  wave. 

Two  .giidl  representations  ar?  useu. 

E(x,z)  = E(^  (x.z)  + i Ej  (x,z)  (Al-a) 

E(x.z)  = * ' * 'x.z)  (A1.b) 

where 

E(x,0)  = E (0)  = 1.0 

The  propagation  LOS  is  assumed  to  sweep  through  the  striated  medium  at  a velocity, 
V,  of  3 X 10'*  cm/sec.  The  field  representation  above  can  be  converted  to  time 
by  setting  x = vt. 

The  field  at  the  receiver  plane  where  z = z^  + z^  is  described  by  a variety 
of  measures:  i;,;  (p^ , (p^*,  x^,  S^,  E^,  Ej,  S^(f),  S^(f),  S^(f),  u^,  a^,  a^,  P(A) 
and  P{((i). 

The  Hg  is  the  e'‘  point  in  the  complex  correlation  function. 

(A2) 

<p^  is  the  variance  of  the  signal  phase.  The  is  defined  as 


w<JO 
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(f)^*  is  the  quantity  used  to  define  mean  square  phase  fluctuation.  Only  in  cer- 
tain situations  is  a real  statistical  quantity.  When  the  mean  square  phase 
fluctuation  is  measurable  and  unambiguous,  it  will  be  described  without  the 
asterisk. 

The  and  are  commonly  used  measures  of  the  amplitude  fluctuations.  The 
is  defined  as 


Both  and  are  determined  from  the  calculated  signals.  Because  of  the  parti- 
cular method  used  to  calculate  x^  and  S^,  the  quoted  values  are  reliable  only  for 
? > 10■^ 

The  and  Ej  are  included  to  provide  an  additional  measure  of  fluctuation 
power.  These  two  measures  are  particularly  useful  for  indicating  gaussian 
quadrature  statistics;  they  also  define  E. 

The  S^(f),  S^(f)  and  Sg(f)  are  the  amplitude,  phase  and  signal  power  spectral 
densities,  respectively.  All  of  these  spectra  are  two-sided  where  frequency 
integrals  are  over  f where  - « < f < co.  The  signal  spectrum  is 

a> 

S^(f)  = y*  [ F“{t)  E (t  + r)  - P]  (A5) 

• 00 

Most  of  these  spectra  will  have  regions  with  a well-defined  power  law  or  gaussian 
behavior  as  they  decrease  at  higher  frequencies.  The  power  law  indices  are 
designated  as  a^,  a^,  and  a^.  The  f-  denotes  a gaussian-1 ike  spectrum. 

The  P(A)  and  P(c())  are  histogram  probability  distributions  for  the  amplitude 
and  phase.  The  Rayleigh  amplitude  distribution  is  plotted  on  each  amplitude  plot 
for  comparison  purposes;  ({>  is  limited  by  + tf  in  calculating  the  phase  distribution. 

Several  of  the  cases  with  large  amplitude  fluctuations  do  not  have  phase 
power  spectral  densities  included.  This  is  because  the  signal  as  shown  in  a 
complex  signal  plot  passes  often  vtiy  rear  the  origin  during  2 deep  fade.  When 
this  happens,  the  phase  measuring  algorithm  has  to  guess  on  which  side  of  the 
origin  the  signal  passed  as  it  processes  the  digitally  represented  signal.  If 
it  guesses  incorrectly,  discontinuities  are  introduced  and  the  spectrum  rapidly 
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approaches  an  f“^  behavior  regardless  of  the  real  spectrum.  The  spectrum  con- 
tains no  useful  information.  If  the  signal  were  analog,  this  numerical  problem 
would  not  exist;  but,  as  shown  earlier  the  signal  passes  near  the  origin  during 
deep  fades  will  still  come  to  dominate  the  real  spectrum. 

This  appendix  is  arranged  in  subsections,  each  containing  the  data  for  a 
particular  assumed  in  situ  environment.  The  subsection  contains  a table  that 
summarizes  the  relevant  parameters  for  the  calculations,  a propagation  space 
plot,  an  Jip/Li  plot,  an  versus  Plot  and  the  plots  unique  to  each  specific 
calculation.  These  sets  are  ordered  in  increasing  An^  for  each  L starting  with 
the  smallest  values  of  L.  The  gaussian  data  are  ordered  by  increasing  instead 
of  An^.  Each  set  of  plots  for  a given  L and  An^  or  z^  has  the  same  figure  number 
indexed  by  letters,  this  includes  a plot  of  the  complex  signal,  amplitude  distri- 
bution, phase  distribution,  amplitude  spectrum,  phase  spectrum  and  complex  signal 
spectrum.  The  phase  spectrum  (letter  designation  "e")  is  ommitted  when  numerically 
ambiguous.  The  last  subsection  contains  the  data  for  the  scattering  layer  thick- 
ness arguments.. 
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Figure  9d.  Amplitude  Power  Spectral  Density 
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Figure  16d.  Amplitude  Power  Spectral  Density 
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Figure  25a.  Signal  Phase  Plot 
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Figure  26c.  Phase  Distribution 
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Figure  27b.  Amplitude  Distribution 
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Figure  28b.  Amplitude  Distribution 
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Figure  30c.  Phase  Distribution 


148 


AFWL-TR-77-183 


*<.08  •!.» 


m iiao  t.0B  z.« 


••M  >^00  tM 

IN  rmsc  cortMENT 


Figure  35a.  Signal  Phase  Plot 
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Figure  35b.  Amplitude  Distribution 
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Figure  35c.  Phase  Distribution 
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Figiire  35d.  Amplitude  Power  Spectral  Density 
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Figure  37a.  Signal  Phase  Plot 
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Figure  38a.  Signal  Phase  Plot 
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Figure  38c.  Phase  Distribution 


193 


96 


AFWL-TR-77-183 


Table  3 
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figure  43a.  Signal  Phase  Plot 
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Figure  44a.  Signal  Phase  Plot 
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Figure  45b.  Amplitude  Distribution 
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Figure  45c.  Phase  Distribution 
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Figure  45e.  Phase  Power  Spectral  Density 
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Figure  46a.  Signal  Phase  Plot 
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Figure  4i)c.  Phase  Distribution 
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Figure  50e.  Phase  Power  Spectral  Density 


243 


PHflSE  COMPONENT 


262 


LOGCfiMPLlTU 
-6.M  -7.20  -•■<»)  -S.eo  -4.80  -4.00 


Figure  53d.  Amplitude  Power  Spectral  Density 
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Figure  54a.  Signal  Phase  Plot 
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Figui^e  55c  Phase  Distribution 
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Figure  61a.  Signal  Phase  Plot 
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Figure  62c.  Phase  Distribution 
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Figure  64b,  Amplitude  Distribution 
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Figure  65a.  Signal  Phase  Plot 
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Figure 


67d.  Signal  Phase  Plot 
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Table  5 

GAUSSIAN  SPECTRUM 
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Figure  83a.  Signal  Phase  Plot 
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Figure  84b.  Amplitude  Distribution 
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Figure  90b.  Amplitude  Distribution 
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Figure  91c.  Phase  Distribution 
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Figure  92c.  Phase  Distribution 
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THICKNESS  DATA 
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Figure  94d.  Amplitude  Dist»'ibution 
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